Multiple light scattering in anisotropic random media 
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In the last decade Diffusing Wave Spectroscopy (DWS) has emerged as a powerful tool to study tur- 
bid media. In this article we develop the formalism to describe light diffusion in general anisotropic 
turbid media. We give explicit formulas to calculate the diffusion tensor and the dynamic absorp- 
tion coefficient, measured in DWS experiments. We apply our theory to uniaxial systems, namely 
nematic liquid crystals, where light is scattered from thermal fluctuations of the local optical axis, 
called director. We perform a detailed analysis of the two essential diffusion constants, parallel and 
perpendicular to the director, in terms of Frank elastic constants, dielectric anisotropy, and applied 
magnetic field. We also point out the relevance of our results to different liquid crystalline systems, 
such as discotic nematics, smectic-A phases, and polymer liquid crystals. Finally, we show that the 
dynamic absorption coefficient is the angular average over the inverse viscosity, which governs the 
dynamics of director ffuctuations. 

PACS numbers: 61.30.-v, 42.70.Df, 78.20.Ci, 78.20.Bh 
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I. INTRODUCTION 

Dynamic light scattering (DLS) is one of our most pow- 
erful probes of the dynamics of materials like sim- 
ple liquids, complex fluids, and liquid crystals. In typical 
experiments, light incident on the sample scatters once, 
and its intensity is measured at a detector. Motion in 
the sample, or more generally fluctuations in the local 
dielectric constant, induce changes in the phase of scat- 
tered light, which give rise to temporal fluctuations of 
the light intensity measured at the detector. Such exper- 
iments probe length scales of order the inverse scattering 
wave vector and hence time scales of order (qv)~^ 
where w is a typical velocity. 

There are many materials such as colloids, emulsions, 
foams, and some liquid crystals that scatter light so 
strongly that the traditional single scattering analysis of 
DLS does not apply. In these materials, light undergoes 
many scattering events before leaving the sample, and 
the transport of light energy is diffusive rather than bal- 
listic. The study of light transport in random or turbid 
media dates back to radiative transfer theory, flrst intro- 
duced as early as 1905 by Schuster These systems are 
characterized by a scattering mean- free-path I, measuring 
the average distance a photon travels before scattering, 
and a transport mean- free-path I* = 1/ {1 — cosds), mea- 
suring the distance beyond which the direction of prop- 
agation is randomized, where 'ds is the scattering angle 
and the angular bracket denotes an average over direc- 
tion weighted by the differential scattering cross section. 
For distances greater than /*, light energy transport is 
described by a diffusion equation with scalar diffusion 
constant D — cl* /3, where c is the speed of light in the 
medium. 

There has been a resurgence of interest in light trans- 
port in turbid and random media Q because of its close 
connection to the problem of Anderson localization ||] 
and electron transport in disordered systems and be- 



cause of the development of Diffusing Wave Spectroscopy 
(DWS) |,|l^,0,[2|jl3|jl|,|l|, which permits useful infor- 
mation to be extracted from dynamic correlations of mul- 
tiply scattered light. Coherent backscattering, a manifes- 
tation of weak as opposed to strong or true Anderson lo- 
calization, has been observed in a number of experiments 
@,0jl|Jl|,|2^,P|2i and discussed in a number of the- 
oretical papers |lff|,||,||,|6|j|^j2|. DWS has opened 
up a whole new field of study. It provides heretofore 
unobtainable information about the the dynamics of tur- 
bid media, including dense colloids [|12|, sheared suspen- 
sions j2^] emulsions [0, and foams ISl]. Because inten- 
sity variations measured at the detector arise from phase 
shifts distributed over many scattering events, DWS de- 
tects dynamic phenomena at much shorter time scales 
than normal DLS. This has permitted the measurement 
of hydrodynamic interaction contributions to the diffu- 
sion constant of colloidal particles and the mea- 
surement of shape fluctuation modes in tense emulsion 
droplets |30j . Photon diffusion and DWS have also found 
application in imaging of objects such a tumors in human 
tissue [|3|. 

With few exceptions |^ , ^ 36,3^], both theory and 
experiment have focussed on diffusive transport and 
DWS in isotropic systems. There are, however, many 
turbid materials such as conventional thermotropic 
and lyotropic liquid crystals, liquid crystalline colloids 
[ p8| , p9[ and emulsions, and also muscle tissues that are 
anisotropic. This paper will develop a general treatment 
of diffusive light transport and DWS in anisotropic media 
with particular applications to nematic liquid crystals. 
Though its inspiration is recent experimental work on co- 
herent backscatter [ ^4p^ ], its purpose is to broaden the 
class of materials to which DWS and its offshoot applica- 
tions like imaging can be applied. A preliminary account 
of this work and experiments on multiple scattering in 
liquid crystals to which it applies were reported in ref- 
erences An alternative derivation of the results 
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reported here and a more detailed account of experiments 
appear in Ref . . A similar treatment of diffusive light 
transport and DWS was developed by Tiggelen, May- 
nard, and Heiderich [^ . 

Anisotropic media differ from isotropic media in two 
important ways: (1) The speed of light depends in gen- 
eral on both the polarization and direction of light prop- 
agation relative to anisotropy axes, and (2) scattering 
cross sections depend not only on the relative direction 
of incoming and scattered light rays but also on their di- 
rection relative to anisotropy axes. Diffusive transport 
in optically active isotropic media with light speeds de- 
pending on the state of circular polarization has been 
studied ||2^. Tiggelen has investigated anisotropic light 
diffusion induced by a magnetic field H through a series 
expansion in H | |3^ . To our knowledge, however, no thor- 
ough treatment of multiple light scattering in optically 
anisotropic media has been published. For electronic sys- 
tems anisotropic diffusion has been studied both theoret- 
ically and experimentally in the context of localization 

Optical anisotropy leads to anisotropic diffusive light 
transport. The equation governing the electric- field au- 
tocorrelation function W{R,T,t) = {E{R,T + t/2) ■ 
eoE{R, T — t/2) ) ( £o denotes the dielectric tensor) is 



d_ 

df 



W{R,T,t) = giR,T) , (1.1) 



where D stands for the anisotropic diffusion tensor. The 
quantity /z(i) is the dynamic absorption coefficient mea- 
sured in DWS experiments ||,|l^,|n|,[2|,[3|jl|Jl| . It re- 
sults from an average of short-time dynamic correlations 
over angle and polarization. We will provide explicit 
formulas for D and ^{t) for general anisotropic systems 
and then concentrate on nematic liquid crystals as one 
example of a uniaxial system. With the preferred axis 
along the unit vector n, the diffusion tensor D reduces 
to £> = D±l + (_D|| — D±)n ® n where 1 is the unit ele- 
ment. The dynamic absorption coefhcient ^,{t) will turn 
out to be the angular average of an inverse viscosity. 

In this paper, we will restrict ourselves to the weak- 
scattering limit, and we will treat multiple scattering via 
the Bethe-Salpeter equation in Sec. III. In isotropic loss- 
less systems, the diffusion equation can be obtained ex- 
actly from the Bethe-Salpeter equation by considering 
only modes associated with the isotropic and "vector" 
spherical harmonics Iqo and Yim- In anisotropic sys- 
tems, all spherical harmonics couple, and the calcula- 
tion of diffusion coefficients involves the inversion of an 
infinite dimensional matrix, which can only be accom- 
plished approximately. We will, therefore, be content 
with a formal expression for the diffusion tensor in gen- 
eral anisotropic media. We will, however, introduce a se- 
quence of approximations to obtain numerical values for 
the diffusion constants in nematic liquid crystals. For- 
tunately, the first term in this sequence undergoes only 
a very small modification in going to the second in this 
sequence. 



Nematic liquid crystals present a difficulty that is 
not generic to anisotropic systems. Light scattering is 
from fluctuations in the direction of the principal axis 
of the dielectric tensor, which is parallel to the local 
Frank director n(r). Fluctuations in n(r) diverge as 
q^^ at small wave number q in the absence of an exter- 
nal aligning magnetic field H . This divergence leads to 
a vanishing scattering mean-free-path for extraordinary- 
to-extraordinary scattering in the limit H ^ 0. The 
diffusion constants are nonetheless well defined and non- 
zero. For, if I tends to zero, scattering takes place almost 
entirely in the forward direction. Thus, light has to un- 
dergo a large number of scattering events before direc- 
tional information is lost, and, as a result, I* is finite. 

The outline of this paper is as follows. In Sec. II, 
we review light propagation including the one-particle 
Green function in homogeneous anisotropic media. In 
Sec. Ill, we treat diffusive transport of light in gen- 
eral anisotropic random media. We introduce the struc- 
ture factor B'^{r, t) to describe fluctuations in the dielec- 
tric tensor, discuss electric-field autocorrelation functions 
and their meaning, and relate them to the averaged two- 
particle Green function. We, then, discuss the one- and 
two-particle Green functions in the weak-scattering limit 
and derive the diffusion equation for light transport from 
the Bethe-Salpeter equation. We introduce the approx- 
imation scheme for the diffusion tensor and look at the 
isotropic limit of our theory. Comments on DWS close 
Sec. III. Sec. IV applies the general results of the pre- 
ceding sections to nematic liquid crystals. A review of 
relevant properties of nematic liquid crystals and light 
propagation in uniaxial media is followed by an explana- 
tion of dielectric tensor fluctuations in nematics. Finally, 
we discuss diffusive light transport and DWS in nematics. 
In particular, we provide explicit numerical calculations 
of the diffusion coefficients D\\ and D±^ as a function of 
Frank clastic constants, dielectric anisotropy, and exter- 
nal magnetic field, and point out their relevance for dif- 
ferent liquid crystalline systems, such as discotic nemat- 
ics, smectic-^ phases, and polymer liquid crystals. The 
numerical calculations are summarized in Figures |7-11. 
They are in excellent agreement with recent experiments 
on the nematic compound 5CB by Jester, Kao, and Yodh 
At the end we address the dynamic absorption 
coefficient. 



II. LIGHT PROPAGATION IN A 
HOMOGENEOUS MEDIUM WITH DIELECTRIC 
ANISOTROPY 

Light propagation in anisotropic dielectric media is 
more complicated than it is in isotropic systems. In par- 
ticular, the electric field is not always transverse, and 
the speed of light depends on polarization and direction 
of propagation. In this section we review light propaga- 
tion in anisotropic media. Following the work of Nelson 
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and Lax [Q, we will introduce sets of polarization vec- 
tors for the electric and dielectric field that will prove to 
be very useful for our forthcoming considerations. We 
start with Maxwell's equations 



and 



div£) 

curli; = - 



1 dB 

c~dt 



divB 







CUrljy = — J ma 

C 



1 dD 



(2.1) 



where Pma and jma are, respectively, the macroscopic 
charge and current densities. We concentrate on a di- 
electric medium with 



D = eoE and B = H 



(2.2) 



and we assume that the dielectric tensor eq is real and 
does not depend on time. Then the energy-balance equa- 
tion reads: 



— U -I- div5 = Jma • E , 

ot 

where we have introduced the energy density 



u = ^{E-£oE + H-H) 

OTT 



(2.3) 



(2.4) 



and the Poynting vector 



Pt{k) = 1 - fe (g) fc 



(2.10) 



is the projection operator on the space perpendicular to 
the unit vector k — k/k. The symbol ® stands for the 
t ens or product: [k ® k]ij = kikj. The solutions to Eq. 
( p.8D provide us with the characteristic light modes of the 
system determined by the refractive index ni{k) — ck/uj 
and the polarization vector ei{k) of the electric field. We 
define the polarization vector of the displacement field by 



d\k) = £oe,{k) 



(2.11) 



For each direction fc, there exist two characteristic light 
modes with associated polarizations d*(fc) _L fc. When 
a plane wave with frequency uj enters the anisotropic 
medium, it splits up into the two characteristic modes 
which travel with different speeds c/ni(k) and wave num- 
bers k = Lojiijc. We also find a third solution with 
Ti'i — CO and e3(fc)||fc, corresponding to a nonpropagating 
mode with w = 0. It violates fe-d'^ = but is necessary to 
construct the complete Green function for wave equation 
( p.6D (see below). 

The two sets of polarization vectors fulfill the biorthog- 
onality condition 



S ^—Ex H 
An 



(2.5) 



Both quantities strongly vary in space and time. If we 
average over one period of oscillation, we obtain their 
averaged values u and S. In the following we use complex 
waves E = E^oe"*"* and H = ffoe"^* whose averaged 
energy densities and Poynting vectors are given by u = 
(Eo ■ eEo + Ho ■ Ho)/{16tt) and S = c{Eo x Ho)/{8n). 

For vanishing sources gima and jma, we obtain the 
homogeneous wave equation for the electric light field 
E{r,ty. 



[curl curl + ^] E{r, t) = 



(2.6) 



Note that all solutions of Eq. (2.(;) have to fulfill 
the trans versality condition for dielectric field waves, 
div£o£! = 0, unless the electric field is static or linear 
in t. Introduction of the plane-wave ansatz 



d\k) ■ e,{k) ^ (z,j = 1,2,3) 



(2.12) 



i. e., they are dual to each other like the basis vectors of 
the real and reciprocal lattice in a crystal. We will use 
both of them as convenient bases for our tensor quan- 
tities throughout the article. To prove condition ( 2.12| ) 
we notice that Pt(k ) and Eq are symmetric tensors and 
derive from Eq. (|2.8|) the condition 



nf nj 



(2.13) 



from which ( p. 12 ) follows after an appropriate normal- 
ization. The vectors a nd d^ are perpendicular to fc. 
Then, again with Eq. (2.8) and the biorthogonality con- 
dition (2.12), one can show that d^ ± d^. In Fig. |^, we 



summarize the geometry for a given propagation direc- 
tion fc. Finally, we recall that in general the refractive 
indices are calculated from Fresnel's equation p7|: 



E{r, t) = Eoe{k) exp[i(fe • r - cut)] 



(2.7) 



in wave equation (2.6) leads to a generalized eigenvalue 
problem 



1 



where 



:Pt(fc)- ^£o]e(fc) =0 , 



1 uj' 



(2.8) 



(2.9) 



3 

E 



Sik; 



(2.14) 



where Si stands for the principal dielectric constants, the 
eigenvalues of Eq, and ki is the component of the unit 
vector fc along the ith eigenvector of £o- This equation 
is equivalent to det[Pt(fc) — So/ri^] = 0. 

Consider now Eq — EQei{k) and Hq, the ampli- 
tude of the magnetic field wave, which are connected via 
Maxwell's relations: 



3 



kx E 







Hq and k X H, 







■£oEo ■ (2.15) 



The second equation also follows from the first once Eq. 
(2^) is solved. For the nonpropagating mode (w = 0), 
Hq, like £^07 is parallel to k. For the propagating modes, 




^fc||e,ffcl 



FIG. 1. The polarization vectors ei{k) and d'(fc) for a given 
propagation direction k in an anisotropic medium. 



Hq is perpendicular to k and So-Eq. Using Eqs. (2.8) and 



15) one shows that electric and magnetic fields carry 
the same amount of field energy {Hq ■ Hq — Eq ■ EqE^) 
and that the averaged energy density iZ of a light wave is 



^0^*0 - ^ 1-^0 1 



(2.16) 



Here the meaning of the scalar Eq becomes clear. It 
is not the magnitude of the electric field, because ei(fc) 
is not a unit vector, instead it basically stands for the 
square root of the energy density of the light mode. The 
poynting vector S does not generally point along k. Its 
projection on k however fulfills the relation 



Q 

k ■ S ~ — u , 

n 



(2.17) 



which is familiar for isotropic systems. 

Finally, we calculate the Gr een function Gq{i — r' ,t — 
t') for the wave equation (2.6), which has a source term 
proportional to the time derivative of jma- In Fourier 
space we have 



[fc2p,(fc)- — £o] Go(fc,c^) = -1 



(2.18) 



We expand Go(A;,a;) in the basis {ei{k)®ej{k)} and cal- 
culate its components by multiplying the last equation. 



respectively, from left and righ t w ith e;(fc) and d''{k). 
Using the eigenvalue equation ( |2.8| ) and the algebra for 
the polarization vectors, described above, we find that 
Go(fc,cj) is diagonal: 



Goik,Lu)^J2 



OJ 



nl{k) 



ea{k) (g) ea{k) 



■ e3(fc) (g) e3(fc) 



(2.19) 



The first term on the right-hand side represents the prop- 
agating part of the Green function. We will always, 
throughout this article, indicate it by Greek indices. It 
acts on the transverse part of the source term. Nelson 
and Lax jisj have calculated an expression for its far field 
in coordinate space, which reduces to a spherical wave in 
an isotropic system. The second term takes into account 
the longitudinal part of the source. In coordinate space, 
it is just a non-propagating delta function. As we go 
on, it will become clear that we can neglect it within 
our approximation outlined in the next section. We will 
therefore skip it when we represent tensors through their 
components. We will also see that for all the tensors 
involved we can neglect the non-diagonal parts. Thus, 
whenever we use a Greek superscript or subscript a it 
refers, respectively, to the basis "vectors" ea{k) (g) ea{k) 
or d" (k) ^ (k) {a — 1, 2). The last quantity we need is 
the inverse Green function [Go(fc, cj)]^^, which wc con- 
veniently represent as 



[Go(fc,^)]-i^^[!^ 



nl{k) 



d°'{k)®d°'{k) 



- — d^{k)(g)d^{k) (2.20) 



by realizing that the unit tensor can be written as 
1 = ^(5;-e,(fe)(g)dJ(fc) =^5^d*(fc)(g)ej(fc) . (2.21) 



III. DIFFUSIVE LIGHT TRANSPORT IN 
ANISOTROPIC RANDOM MEDIA 

We now consider a dielectric medium in which the 
dielectric tensor possesses a randomly fluctuating part 
5e{r,t) in addition to the homogeneous term eo- We 
treat <5e(r,i) as a Gaussian random variable with vari- 
ance characterized by 



\iN) 



(3.1) 



where uj is the frequency of light. The superscript (N) 
means that we interchange the second and third index 
in the tensor product Se (g Se to define -B": [B'^]ijki oc 
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(feifc Ssji ). We call B'^{r, t) the structure factor of the 
system. It is measured in single light scattering experi- 
ments [|] and contains information about the elastic and 
dynamic properties of a system. 

The electric light field E{r,t) obeys the homogeneous 
wave equation 



1 92 

[curl curl + (£„ + Se{r, t)) — —] E{r, t) 







(3.2) 



where we have used an adiabatic approximation to pull 
i5e(r,t) in front of the time derivatives. It is valid if 
i5e(r,t) varies on time scales much longer than the time 
period of light and the passage time of a light ray through 
an inhomogeneous medium. 

O ur t ask is to calculate measurable quantities from 
Eq. (3.2). A very general one is given by the spatial and 
temporal autocorrelation function for the electric light 
field, {E{ri,ti)^E{r2,t2)), which we write with the help 
of center-of- "mass" {R, T) and relative (r, t) coordinates: 



W{R,r,T, t) 



E 



T 



(3.3) 



It is a second rank tensor. For r = we derive the scalar 

W{R,T,t)^(^E(^R,T+^-)-eoE*(^R,T-^-)) , 

(3.4) 

which at t = is equal to the energy density of the 
light field when the small fluctuating part is neglected. 
The T dependence is, e.g., due to time-modulated light 
sources, and R describes variations of the energy density 
on long length scales. For t 0, W{R,T,t) reflects the 
dynamics of the scattering medium through its depen- 
dence on the structure factor B'^{r,t). It is measured 
in dynamic light scattering experiments either in single 
scattering g or with DWS Wc wiU see below that 
the Fourier transform with respect to r, W{R, k,T,t), 
is diagonal in the basis {ea{k) (g) e^(fc)} on length {R) 
and time (T) scales much longer than the wavelength and 
time period of light, 



W{R,k,T,t) 



'W°'{R,k,T,t) e^k) (g) ea{k) 



(3.5) 



For t = we can, therefore, interpret W°'{R,k,T,t) as 
the energy density of a light wave with wave vector k and 
polarization ea{k). To clarify this interpretation we look 
at the Fourier transform 



/ / K n\ ,/ K n 

E(k + -,. + -)^E*(k--,u~- 



(3.6) 



Since R and T describe variations on long scales we 
have K <^ k and ft uj. Thus the amplitudes 
E{kzLK/2,ujzLn/2) are srongly peaked around u and k 
which can be identified, respectively, with the frequency 
and wave vector of the light waves in the medium. 

To calculate the autocorrelation function W{R, r, T, t) 
for special light sources Wo{R,r,T,t) and/or given 
boundary conditions, we need the averaged "two- 
particle" Green function 



(3.7) 



where and denote the retarded and advanced 
Green functions with G'^ = [G^] * . Then 



Wil) = y"d2$(l,2)W^o(2) 



(3.8) 



(for notation see Appendix A). In the second subsec- 
tion we will derive the diffusion pole of i.e., that 
part of $ which corresponds to a diffusion equation 
for W{R, k,T,t) in the variables R and T. It governs 
the propagation of light energy at long length and time 
scales. The time t will appear in an absorption term of 
the diffusion equation, which is zero for t = 0. In the 
first subsection we calculate the averaged retarded and 
advanced one-particle Green functions, which we need for 
the derivation of 



A. The averaged one-particle Green function 

The one-particle Green functions (G^) and (G^) fol- 
low frompDyson's equation | |50| ], which we give in a formal 
notationo 



(G^/^) = Go + Go (G^/-^) 



(3.9) 



We have introduced the retarded and advanced second- 
rank-tensor mass operators and S^. In the weak- 
scattering approximation, where the elements of the 
structure factor B'^{r,t = 0) are assumed to be much 
smaller than one, these mass operators are proportional 
to B'^{r, t = 0) and can be written in frequency space as 



S^/^(r,t^) = B'^(r,t 0){G^/^){r,u;) 



(3.10) 



Figure |2| gives a diagrammatic representation of 
S"^/" ^(r, uj). Finally, in momentum space, we get from 
Eq. (pj), 



*In the following we will use a coordinate-free representa- 
tion for tensors and their contractions which we explain here. 
If A and B are second rank tensors so is AB with compo- 
nents [yl-B]ij — '^1^ AikBkj ■ For fourth rank tensors C and 
D we form a tensor CD of the same rank with components 
[CD]ijki = J2mn CijmnD^nki- Finally, CA is a second rank 
tensor and has components [CAJi^ = ^^.^ CijkiAki- 
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(3.11) valid in the limit 0. We finally obtain 



The inverse Green function Gg ^{k,uj) is diagonal in the 
basis {rf*(fe) (g) d^(fc)}, and S^/"^(fc,w) causes a small 
pertubation. From perturbation theory we know that to 
zeroth order in S^/"* the eigenvectors of GQ^{k,uj) — 
I]^/"^(fc,a;) are unchanged, whereas to first order in 



2 



(3.17) 



where we sum over all possible scattering events of the in- 
coming light mode {fc" | ea{k)} into modes {q^ \ ep{k)}. 
The single scattering event is described by the structure 
factor 



-<5e(O,0) 



^0) 



(5e(r,0) 



FIG. 2. The mass operators S^''^(r,a') in the weak-scat- 
tering approximation consist of two scattering events from 
inhomogeneities in the dielectric tensor, which are tied to- 
gether by spatial correlations [_B"(r,t = 0)] and the averaged 
propagators {G^^'^){r,uj) of the electric field. 

i^j-^g diagonal elements 

[S^/^(fc,c^)],; = e,(fe) • S«/^(fc,c^)e,(fc) (3.12) 

contribute to the eigenvalues of {G^^^) renormalizing the 
wave numbers ^ n, of the light modes. The significant 
effect of H^^^ comes from its imaginary part. Then, in 
its final form the propagating part of the Green functions 
reads 



(G«/^)(fc,c.)f 
with 

[{G^'^){k,Lo)r 



2 

El 



;G^/^)(fc,c.)]"e,(fc)®e„(fe) 

(3.13) 



nl{k) cna{k)la{k,uj) 



(3.14) 



We have introduced the scattering mean-free-path 
la{k, oj) of the light mode {A;" | ea{k)} which travels with 
a wave vector fc" = —Uak and polarization ea{k): 



la{k, Uj) 



■rv(fc)[ImS^(^n„fc,w)]„ 



(3.15) 



To arrive at la{k,uj), we used a CPA approximation re- 
placing the argument fe of [S-'*-/'^ (fc, uj)\a by the wave vec- 
tor ^ riak of the li ght m odes where Gq diverges. ImS-'^ 
follows from Eq. ( |3.10 ) using the momentum shell ap- 
proximation 



[Im(G«/^)(fc,i 



-i,-n^{k)5{^n^{k)-k) (3.16) 



\B 



<{t = 0)]ai3 ^ ( \ea{k) ■ 6£{qs,0)ep{q) 



(3.18) 



which is proportional to the differential scattering cross 
section. We also introduced the scattering vector 



{riak - riisq) 



(3.19) 



In Eq. ( 3.17 ) and throughout the article, we use a short- 
hand notation for angular integration: 



(Kin 



(2^) 



(3.20) 



We point out that the scattering mean-free-path la{k, oj) 
now depends on the direction k and the polarization 
ea{k) of the light mode. Following the work of Nel- 
son and Lax p8[ , it is straightforward to show that 
la{k,Lij) determines an exponential decay of the far field 
of {G^){R,uj) travelling in the direction k with polar- 
ization ea{k). In our analysis, we have neglected the 
off-diagonal components of S^/'^, so that the polariza- 
tion vectors of our light rays are identical to those of the 
homogeneous medium. We will use this approximation 
in the next subsection. We conclude with the definition 
of the two quantities 



(G^Xfc^ 



(G^)(fc_,c._) 



AGt{K,n) 
and 

AS^(K,f]) = S^(fe+,w+) - S^(fe_,w_) 
where 

k±=k±K/2 and lu±=uj±Q/2. 



(3.21) 
(3.22) 
(3.23) 



We will use both of these in the following, especially when 
K = and = 0: 



AG^(O,0) = 2fIm(G^)(fe,t^) 



(3.24) 



and 



AS^(O,0) = 2iImS^(fc,tj) 



2i 



c na{k)la{k,uj) 



(3.25) 
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B. The averaged two-particle Green function 

The averaged two-particle Green function $ obeys the 
Bethe-Salpeter equation ||50t| : 

(3.26) 

where U stands for the irreducible vertex function or 
intensity operator, which in the weak-scattering approxi- 
mation, equals the structure factor: U k, B. The Bethe- 
Salpeter equation is best handled in momentum and fre- 
quency space. With all arguments, the Green function 
is rj,t) (see Appendix A). The variables K, Q. 

correspond to the center of "mass" coordinates i2, T, 
introduced in the introduction to this section, and the 
wave vectors fc, fe' to relative coordinates r, r'. The su- 
perscript uj is the light frequency, and the t dependence 
exphcitly comes from the structure factor B^j,,(t). The 
explicit form of the Bethe-Salpeter equation reads 

= mK,nnit , (3.27) 

which we derive in Appendix A. In the last equation we 
have introduced an abbreviation for the tensor product 
of the averaged one-particle Green functions: 



ft{K,n) = [{G'''){k+,iU+)® (G^)(fc-,c._)]' 



(3.28) 



and a shorthand notation 1^*^^, for the combination of 
the delta function and the unit element of fourth rank 
tensors: 

[^itkki (27r)3 d{k ~ k') {S.kSji + SuS,k)/2 . (3.29) 

Here i,j,k and I are Cartesian indices. The Bcthe- 
Salpeter equation ( [3.27 ) can be solved formally by it- 
eration, leading to a sum of ladder diagrams (see Fig. 



k' 



k' 



k' 



k. 



+ 



k' 



+ 



FIG. 3. The two-particle Green function $ as a sum of lad- 
der diagrams. f%{K,Q.) propagates two electric field modes 
with wave vectors k'^ and k'_. Both are scattered by the 
same inhomogeneity in the dielectric tensor described by the 
structure factor B'^y{t) for the scattering vector k ~ k' . 



The multiple integrals and the sum can be done ana- 
lytically for (5-function correlations for — > 0, and 
the solution represents the diffusion pole of For the 
anisotropic, long-range correlations of our problem this 
procedure is not possible. VoUhardt and Wdlfle jsi] de- 
rived the diffusion pole for isotropic electron transport 
directly from Eq. (3.27), and MacKintosh and John 
applied their method to light. In the anisotropic case one 
has to be more careful, as Wolfle and Bhatt showed for 
electrons Q]. We will proceed a similar way, however 
our problem is more difficult because we have to deal 
with the different polarizations of light. 

Let •^^^\K,n,t) and X('^\K,n,t) be, respectively, 
the nth eigenvector and eigenvalue of the integral opera- 
tor of Eq. (|3^), 



-/^(K,17)B-,^(0]*i:^ = At")*!") , 

(3.30) 

and '^j. {K, f^, i) the eigenvectors of the Hermitian ad- 
joint operator 153]. Then it is straightforward to show 
|l5| , with the aid of the completeness relation 



(3.31) 



that 



■rk'iK,n) (3.32) 



solves the Bethe-Salpeter equation. When K = and 
n = t = 0, the quantity AGy (O,0) , defined in Eq. 
( |3.2l|) , is an eigenvector of Eq. ( 3.30 ) with eigenvalue 
-V^(0,0, 0) = 0. To prove this statement we use a spe- 
cial case of one of the Ward identities (see Appendix |B l| ), 

AS^(O,0)= J j^Bt,,{t^0)AGUO,0) , (3.33) 



and the relation 



ag;(^(o,o) = /^(o,o)as-(o,o) 



(3.34) 



The first equation is obvious from the definition of the 
mass operator, and the second one is given in Appendix 
B2. We have identified the diffusion pole, as we shall 



explicitly see soon. The result is valid beyond the weak- 
scattering approximation and based o n the Ward identi- 
ties All other eigenvalues of Eq. ( 3.30| ) are positive, 
and in the real-space coordinate i?, they give exponen- 
tially decaying contributions to ^f,i.,{K, ft, t), which are 
not important at long length scales [|l5| . 

The diffusion approximation follows when we calculate 
;\(o) (^j^^ with the help of perturbation theory in the 
limit K,n,t 0. First, we have to get the eigenvalue 
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equation in this limit. To proceed we n eed the following 
equation, which we derive in Appendix |B 2|: 



dk 



K 



2uj 



(3.35) 



As usual, the Greek superscript or subscript a refers, 
respectively, to the basis "vectors" ea(fc) <Ei ea{k) or 
d°'{k) (E) d'^iit). Equation (|!3|) gives [AG%{K,n)]°' 
to first order in K, fl and S. It shows clearly where 
the wave vector K and the frequency Q come in. Now, 
we formulate eigenvalue equation ( 3.30| ) in components, 
choosing K — and = 0. Then, we apply the last 
equation and arrive at 



Asno, o)l„ - sail . jf - 1; n I *; 



( 



- [AG^(O,0)]"^ 



(2^ 



mkAt)UK (3-36) 



as the eigenvalue equation, which is correct up to first 
order in K , Q, and S. In a next step we use the ansatz 

vI/^cx[AG^(O,0)]"[*ovr + ^*.<(fc)] , (3.37) 



for the eigenfunctions to turn the eigenvalue equation 
into a matrix equation. The first factor on the right- 
hand side forces k to equal —na{k) in the momentum 
shell approximation: 

[AG't(O,0)r K^in-nJk)6(^nJk)-k) . (3.38) 

The amplitude V&o represents the zeroth eigenvector, and 
the second term involves a complete set of real basis func- 
tions ipi{k) on the unit sphere and for the space spanned 
by ei (g) ei and 62 (X" 62. We will comment below on how 
to choose these basis functions. Keeping only dominant 
terms in the small fl, K, and t limit, we can write the 
eigenvalue equation in matrix form as 



-iQ{K) 



B 



So Si 



(3.39) 



The matrix B basically represents the structure factor 
(0)]q/3 in our chosen basis. 



\B] 



a, 13 
-if. 



k"(fcX(fe) 



r(fc)^f(9)] [Bt.,,{Q)U} 



(3.40) 



The vector Q with components 



^ -Ik" 



dk 



Kftik) (3.41) 



is linear in K. Another important term is the coefficient 

t3 



(3.42) 



a,/3 ■ 



where 



— 1 



= ^E/^"'^'^'(^) • (3.43) 

a 

For reasons that will become clear shortly, we will refer 
to ^{oj,t) as a dynamic absorption coefficient. It is an 
average over dynamical modes that decay in time. As 
a result, it increases from zero with increasing t > 0. 
Finally we have the constant 



E 



d^k 



(3.44) 



(27r)3 



[AS^(O,0)]„[AG^(O,0)]" . (3.45) 



The second equation above states that S is propor- 
tional to the normalization factor of the eigenfunction 
AG^(O,0), with AS^(O,0) being the eigenfunction of 
the Hermitian conjugated problem. The vector Sq and 
the matrix Si are irrelevant here. For K,H,t = the 
zeroth eigenvalue A*^*^-* corresponds to th e eigenvector 
[1, 0, 0, . . .] in the eigenvalue equation (3.39), and we want 
to know how it evolves for small K, and t. The vector 
G{K) couples the 00-element with the matrix B. The 
coupling can be removed by an orthogonal transforma- 
tion with 



1 i[giK)YB-^ 
-iB-^g{K) 1 



(3.46) 



which essentially renormalizes the 00-element. Calculat- 
ing A*-"^ from the transformed system to leading order 
in shows that A'-^-'S is equal to the renormalized 00- 
element, or 



A(") = ^ 



-in + ^l{uJ, t) + ^ g{K) ■ B-^g{K) 



(3.47) 



Now, we are ready to write down the diffusion approxi- 
mation of the averaged two-particle Green function from 
Eq. (pi. 



(Knt)~ 1 AG^(O,O)0AG^,(O,O) 
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with 



K-D{uj)K = ^g{K)-B-^g{K) (3.49) 



and N = -2n^u;'^ /{nc^). When fj.{uj,t) = the Green 
function ^'^i.,{K ,il,t) possesses a simple diffusion pole 
with an anisotropic diffusion tensor D{lo). In addition 
we have the dynamic absorption coefficient ^{oj, t), which 
forms the basis of DWS. We will discuss it in the next 

subsection. 

From Eq. ( 3.49| ), it is clear that the problem of calcu- 
lating D{u}) reduces to the problems of calculating Q{K) 
and the inverse matrix . If we had at our disposal 
functions ipi(k), that diagonalize B, our task would be 
simple: 



K ■ D(uj)K ^ V 



[G{K)]t 
[Bh 



(3.50) 



Unfortunately these functions are difficult to find in 
anisotropic systems, and we are reduced to seeking us- 
able approximation schemes. One scheme is to choose 
functions <^j(fc) so that only a few components of G{K), 
say [Q{K)]i for i — 1,2, are non-zero. Then, if the off- 
diagonal components of [B]ij, coupling j for j > 2 to 
i ioT i = 1,2, are small, [g{K)]i[B~^]ij[g{K)]j can be 
approximated by using only the i,j ~ 1,2 submatrix of 
\B\ij to calculate \B^^\ij. We will use this approach in 
numerical calculations of D{lo) in nematic liquid crystals 
in the next section. 

To test our theory we apply it to isotropic systems. 
We concentrate on the case, where fluctuations in the 
dielectric tensor originate from density fluctuations only: 
5e (X 5q\. Then the components of the structure factor 
assume the form 

[S^»g^(0)]a/3 = [ea(fe) • ep{q)f i?-(cosi?,) , (3.51) 

where B'^{cos'ds) solely depends on the scattering an- 
gle -ds- In isotropic systems, all the polarization vectors 
perpendicular to k are equivalent, and the diffusion con- 
stant should not depend on the special choice we make. 
Therefore, to treat the factor [ea{k) ■ e^(q)]^, we can add 
integrations / (i0i/(27r) J (i02/(27r) in the formula for B. 
The angles (pi and (/)2, respectively, describe rotations of 
the polarization vectors about k and q. The integration 
over 01 and 4)2 for a fixed scattering angle "dg is straight- 
forward, and we obtain 



27r 



27r 



[ea(fc,0i) ■ ep{q,cj)2) 



1 + cos^-ds 
4e§ 



(3.52) 



where the dielectric constant Eq comes in through the 
normalization of the polarization vectors. With this 



trick we are able to show that E is essentially di- 
agonalized by spherical harmonicsl. As basis func- 
tions we choose ipiim{k) = [y;™(t9, <^), yim(z9, <^)] and 

f2hn{k) = [yim(t9, <^),-F/m(l9,(^)] with C0S1? = fc • 

K. We also use the addition theorem, P;(cosz?s) oc 
E„^YLi^,mm{^',^') and arrive at [Br^;\? - 



The only non-zero clement of Q is [tJJno = 
— 2£oif/\/3, and we finally find the formula for the diffu- 
sion constant: 

U = TTCn 

3 " 



dcosdsil - cosi?s)(l + cos^ds)B'^{cosds) 



(3.53) 



where cq stands for the speed of light in the system. 
Compared to the scalar case we have an additional fac- 
tor 1 -I- cos^-iJs- B'^{cos'ds) alone describes a scattering 
process with incoming and outgoing polarization perpen- 
dicular to the scattering plane (VV scattering), whereas 
cos^dgB'^icofids) belongs to scattering with polarizations 
in the plane (HH scattering) In a group theo- 

retical language [B^^^f,{Q)\ap transforms under a high- 
dimensional identity representation of SO{2>). The rela- 
tion cx Sii'Smm' then means that we have de- 
composed [S]y™„j/ via the irreducible representations of 
S0{3). For less symmetric phases, like a nematic liq- 
uid crystal, we can at least partially diagonalize B with 
the help of group theory. The relevant symmetry group 
Daoh has only one-dimensional representions induced by 
the functions exp{imip), where if is the azimuthal angle 
around the symmetry axis. Decomposing B with the help 



f'l'm' 



where dif- 



of spherical harmonics gives [B 

ferent I and now couple to each other (see nex t sec tion) . 
The numerator of the Green function in Eq. (^7 



tains an interesting effect. The second factor AG^, (0, 0) 
collects the light sources or the incoming light waves. The 
ffi-st factor determines the energy density W°'{R, k, T, 0) 
of an outgoing light wave independent of the sources. 
This means that in the diffusion approximation the 
outgoing light looses all its correlations with the light 
sources. Integrating [AG^, (0, 0)]ct over the wavenum- 
ber k shows that the ratio for the energy densi- 
ties W°'{R,k,T,0) in polarization states 1 and 2 is 
[ni(fc)/n 2(fc)] '^. Experiments measure light intensities. 
In Eq. (|2.17D we learned that only the projection of 
the Poynting vector S" on fe is simple: S" ■ k = 
{c/na)W°'{R,k,T,0). Its ratio for the output polariza- 
tions 1 and 2 is [ni(fc)/n2(fe)]^ which gives the ratio of the 
output intensities when the Poynting vectors are parallel 
to fe. This is the case for light travelling along the prin- 



^Our formulas are written for real basis functions, a gener- 
alization to complex ones is straightforward. 



9 



cipal a xes of the dielectric tensor. The Green function 
in Eq. ( 3.48 ) suggests that there is a diffusion equation 



t)-V ■ D{lu)V] W{R, . , t) = q{R, 



(3.56) 



for each hght wave with direction k and polarization a. 
This does not mean that we have additional conserved 
quantities besides the energy density. It only means that 
after randomizing the incoming light the distribution of 
the light modes in the light field stays the same. 



C. Diffusing Wave Spectroscopy 

If we sum over the two polarization states and int egrate 
over all wave vectors k the Green function in Eq. ( 3.4§| ) 
is equivalent to a diffusion equation for the scalar time- 
correlation function W[R,T,t), introduced in Eq. (^.4[), 



d_ 
df 



V ■ D{uj)V + ^l{LO,t) W{R,T,t) ^ g{R,T) 



(3.54) 



This equation is the basis of Diffusing Wave Spectroscopy 
(DWS) [|l^. Solving it for special boundary conditions 
and sources g{R,T), which depend on experimental ar- 
rangements, gives W{R,T,t) in terms of the dynamic 
absorption coefficient fi{uj,t). In our derivation of the 
diffusion pole we had to restrict ourselves to times t with 
[B-„^,(0) - B",^^,{t)U « [B",^^AO)U. In this time 
range we expect being able to perform a Taylor expansion 
in t which gives a linear time dependence for fi(uj,t): 



(3.55) 



The constant fj,o reflects some averaged dynamical prop- 
erties of the system. For the diffusion of particles in 
colloidal suspensions, fig = 2Dbi^^/{1*c), where Db is 
the self-diffusion constant of the particles |4^ . The con- 
dition just imposed on t means that DWS probes the 
dynamics of a system on a much shorter time scale than 
single light scattering does, which requires times where 
[B'^^^0{t)]af3 has already decayed considerably. DWS, 
therefore, offers the possibility of studying the validity 
of the dynamical description of the system on short time 
scales. This was done for colloidal suspensions by Kao et 
al. who studied the Brownian motion of single par- 
ticles at short times where the mean-square diplacement 
is not simply proportional to t p9| . As a resu lt, fi{u!,t) 



does not follow the linear time law of Eq. ( 3.55 ). Experi- 
mentalists prefer a different picture for DWS which they 
have developed for colloidal suspensions |^^. They sum 
up all possible light paths in the scattering medium to 
arrive at the time-correlation function after some aver- 
aging procedure. We will show here that our approach 
is totally equivalent to this picture. However it has the 
advantage that it automatically tells us how to perform 
this averaging procedure. 

DWS experiments are usually performed with contin- 
uous light sources and the diffusion equation reduces to 



We can rewrite this equation as the Laplace transform of 
a problem, where the source is a light pulse. 



r_9_ 

97 



-V ■D{lo)V P{R,T)^d{T)g{R, .) (3.57) 



with 

/oo 
P{R,T)exp[-fi{uj,t)T]dT . (3.58) 
-0 

The lower limit —0 means a small negative time r in order 
to pick up the delta function. We can interpret P{R, t) 
(after an appropriate normalization) as the probability 
that light, emitted by the source at time t — Q, arrives 
at the detector at point R after a time r. Then the time- 
correlation function W{R, . , t) follows after a summation 
over all light paths where each light path contributes a 
factor exp[— /.t(u;, i)T] to the decay of W{R, .,t). For 
isotropic systems r is directly connected to the path 
length s = cr/n of light and P{R,t) also represents the 
path-length distribution. The exact form of P{R, t) de- 
pends on the choice of the light source and the boundary 
conditions |49| which we don't address here. 



IV. LIGHT DIFFUSION IN NEMATIC LIQUID 
CRYSTALS 

The nematic liquid crystalline phase consists of rodlike 
organic molecules which tend to align parallel to each 
other but which show no long-range positional order of 
their centers of mass. The local average direction of the 
molecules is described by a unit vector n{r,t) called di- 
rector. It appears in the local dielectric tensor 



£(r, t) =e±l + Ae[n{r, t) ® nir, t)] 



(4.1) 



where £j_ and ey are the dielectric constants for electric 
fields, respectively, perpendicular and parallel to the di- 
rector and where Ae = ey — £±_ stands for the dielectric 
anisotropy. The energetically favored state of a nematic 
phase is a uniform director field n{r,t) = rip through- 
out the sample. Its distortion costs energy which can be 
calculated from the Frank-Oseen-Zocher free energy [Q: 



- j [ifi(divn)2 + K2{n ■ cnvXnf 
Jr Ksin X curln)2 - Ax(n • H)^] (f 



'r . 



(4.2) 



There are three characteristic distortions, called splay 
(Ki), twist {K2) and bend {K3), where Ki, K2 and 
K3 are the Frank elastic constants. We also include a 
magnetic-field term with = X|| ~ X-L the anisotropy 
of the magnetic susceptibility. If A^ > an alignment of 
the director parallel to the field H is favored. Even in a 
uniformly aligned sample there exist thermally induced 
fluctuations of the director: 
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n(r, t) = no + 5n{r, t) 



(4.3) 



a 



Ae/e_L 



(4.8) 



They lead to fluctuations in the local dielectric tensor and 
hence scatter light. This is the physical phenomenon for 
which we want to calculate the diffusion approximation 
of light. In order to apply the formulas of the preceed- 
ing section, we have to look first at the light propaga- 
tion in an uniaxial media. Then we need to calculate 
the structure factor, associated with the director fluctu- 
ations, which governs the single- light-scattering event. 



A. Light propagation in uniaxial media 

Let us, for a moment, suppress the fluctuations of the 
director and look at a homogeneous medium with a uni- 
axial dielectric tensor 



£o = £±1 + Ae[no ® no] 



(4.4) 



The equilibrium value Uq of the director is also called the 
optical axis because it establishes a special axis for light 
propagation. The two light modes follow from solving 
the eigenvalue equation ( |2.8| ). They are well described 
in the literature We mainly summarize the results 
here, introduce some notation, and perform an interest- 
ing variable transformation for later use. 

One light mode is immediately obvious. Its polariza- 
tion vectors 62 (fc) and d?{k) are both perpendicular to 
no and k with a refractive index 712 = -v/^T- Since it be- 
haves as in an isotropic system, it is named the ordinary 
light mode. We choose no parallel to the z axis and write 
k in spherical coordinates: 



no 



and k 



(4.5) 



Then the ordinary light ray is represented by 



"2 = ^nd 62 (fc) = — 

"2 



-smi^fc 
cos^Sfe 




(4.6) 



In the extraordinary light mode, the refractive in- 
dex depends on k, and the polarization vector ei(fc) 
is generally not perpendicular to fc. It follows from 
ei{k) = eQ"^d^(fc) where the polarization vector d}{k) 
is determined by d}{k) _L dP{k) and d^{k) _L k. The 
refractive index can be calculated from the eigenvalue 
equation (p^) or Fresnel's equation (2.14): 



ni{k) — n2ni{k) with ni(fc) = 



1 



-,1/2 



1 + acos'^-dk 



(4.7) 



where we have introduced the relative dielectric 
anisotropy of the system. 



Equation (4/7) does not represent the most symmetric 
form for ni{k), but we found it useful for our calculations. 
Finally, we get 



ei(fc) 



ni{k) 



"2 



COSifk 

sin.^fc 




sini?i 



(4.9) 



In the structure factor, k always appears together with 
the refraction index. For the extraordinary mode we, 
therefore, replace the angular variables cos'&k and sin'&k 
by an equivalent set, 

Ck ■= ni{k) cos-dk and 5*^ := ni(fe) sirnJ^, , (4.10) 
in which the refraction index takes the form 



ni(fe) = ^1 + ail -CD 
The "trigonometric" identity, 

Sl + {1 + a)Cl = 1 + a 



(4.11) 



(4.12) 



is valid. The new coordinate Ck ranges from —1 to 1 and 
contains the same information as cosdk- Thus, we can 
choose, e.g., spherical harmonics as basis functions on the 
unit sphere with cou'dk replaced by Ck- The usefulness of 
Ck appears when we calculate its differential with respect 
to cos-dk'- 



(1 -f a) dCk = n\{k)dcos-dk 



(4.13) 



The differential dCk absorbs the factor n\{k) which al- 
ways comes with our integrals thus making them easier: 



nl{l + a) 



dCkdipk 
(27r)3 



(4.14) 



Alternatively, we can say that we have constructed a com- 
plete set of basis functions in the variables cos'dk and ipk 
with respect to the weight function ni(fc). In the fol- 
lowing we will also use the notation Ck — cosi^fc and 
Sk = simJfc for ordinary light modes. 



B. Structure factor for director fluctuations 

Here we will derive the structure factor for director 
fluctuations. To begin we recall the general form of the 
time-dependent structure factor for light scattering in an 
anisotropic dielectric. 



'(*)]a/3 



UJ 



[ea{k) ■ Se{qs,t)ep{q) 
[eaik) ■ 6e{qs,0)ef3{q)Y 



(4.15) 
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where Qs = ^ {uak — n/sq) stands for the scattering vec- 
tor. We first need the Fourier component 6£{qs,t) of 
the fluctuating part of the dielectric tensor. We insert 
n(r, t) — no + Sn{r, t) into e(r, t) of Eq. (4.1) and collect 
the first-order terms in Sn{r,t). After a Fourier trans- 
formation we arrive at 

S£{qs,t) = Ae [no ® 5n{qs,t) + 6n{qs,t) ® no] (4.16) 

with 6n{qs,t) the amplitude of a director mode. It has 
only two components since the director is a unit vector. 
Furthermore, for small fluctuations, 6n is perpendicular 
to no. An appropriate basis for 6n{qs,t), as shown in 
Fig. 1^, consists of a unit vector Ui{qs), lying in the 
plane defined by no and qs, and a second one, U2{qs), 
perpendicular to this plane: 

Sn{qs,t) ^ STii{qs,t)ui{qs) + Sn2{qs,t)u2{qs) ■ (4.17) 



*■ U2{qs 



Ui{qs) 

FIG. 4. Basis vectors ui{qs) and U2{qs) for given director 
no and wave vector q^. 

Next we need the temporal autocorrelation function of 
the director modes [Q, which is rather complex. We 
state it and then explain the individual terms: 



{dniqs,t)(g,Sn*{qs,0)] 



<5=1 

with 



2 

^ Ksiqs) 



kBT 



exp 



Ksjqs) 



usiqs) (gi usiqs) (4.18) 



Ksiqs) 



Ksql 



(4.19) 



where q± and are the components of q^, perpendic- 
ular and parallel to no. The chosen basis for Sn{qs,t) 
provides the normal coordinates, because the correlation 
function is already diagonal. The free energy of one direc- 
tor mode follows from the Frank-Oseen-Zocher free en- 
ergy dU) as Ks{qs)\Sns{qs,t)\'^. For a general vector qs 
it is either composed of splay and bend (5 = 1) or twist 
and bend distortions (i5 = 2). The factor k-^T / Ks{qs) 
results from the equipartition theorem giving the mean- 
square amplitude of each mode. The dynamics of the 
director modes is described by the Leslie-Erickson equa- 
tions Ba]. They combine the Navier-Stokes equation for 



a uniaxial media with dynamical equations for the direc- 
tor. A detailed analysis |^ shows that director modes 
are diffusive with a relaxation frequency given by the quo- 
tient of elastic [Ks{qs)] and viscous [775(^5)] forces. This 
is the origin of the exponential factor in the correlation 
function. The viscosity rjs{qs) is a combination of sev- 
eral Leslie viscosities which we will address in subsection 



IV D. We neglect a fast relaxing part which comes from 
the coupling between the director and the fluid motion. 
Finally we are able to write down the structure factor 



EN{ea,ei3,us) 
TT-TTS exp 



■5=1 



Ksiqs 



Ksjqs) 
VsiQs) 



(4.20) 



with 



N{ea, Bp, Us) ^ [(no • e^)(M5 • e^) -f {iis ■ ep){nQ ■ e^)]^ 

(4.21) 

a geometrical factor. It has two interesting implications. 
First, there exists no scattering of an ordinary light ray 
into an ordinary light ray because 62 -L no, and therefore 
= 0. Second, forward and backward scattering along 
no is always forbidden. The other terms of [^^^^^(t)]^^ 
are familiar from the previous discussion. The scatter- 
ing mean- free-path la{k,uj), defined in Eq. ( 3.17 ), has 



been already discussed in detail by two groups [ l56| , p7t . 
Its dependence on k and polarization was calculated by 
Val'kov and Romanov [^ . We just stress one point. The 
structure factor diverges for iJ = and — > 0, and we 
expect the scattering mean-free-paths to be zero in an in- 
finite medium. But = can only occur in scattering 
events where the extraordinary polarization is preserved. 
Hence, only the scattering mean-free-path li{k,oj) of an 
extraordinary light-ray is zero. In Appendix C we give 
[B^i^i(0)]ii for q^ k^ (the notation we use there 
is explained below). From this form it is obvious that 
li^{k,uj) diverges weakly like —\tl[1^xH'^ / {K^nluj'^ / c^)] 
for magnetic fields much smaller than K3 / Axn2'-^ / c 
The explicit form of the structure factor is very complex 
because it explicitly depends on the direction of the in- 
coming (fc) and outgoing (q) light through the geometri- 
cal factor. We mention some symmetry properties which 
facilitate its handling. The first two are quite obvious: 



[Bt.^,{0)U^[B-,k^{0)h 



and 



[Bt.^,{0)U^[B-^,..,,{0)U ■ 



(4.22) 



(4.23) 



The structure factor has to refiect the symmetry of the 
nematic phase described by the group Daoh- For exam- 
ple, it has to exhibit the rotational symmetry about the 
director no- If we choose 
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siwdk cosifik 
sin-dk simpk 

COS'&k 



and q 



, (4.24) 



then (0)]a,3 depends only on the relative azimuthal 

angle ip — ipq ~ ipk- The existence of mirror planes, con- 
taining Uq, implies that the structure factor should be 
invariant under (p — > —ip. In appendix C we give its ex- 
plicit form in terms of Ck, Sk, Cq, Sq and ip. We also 
introduce a scaled structure factor [Bf.o,qi3{0)]ai3 through 



[Bt.q,{0)U ^ (As)' -J 



kBT 



[Bk^q^U , (4.25) 



which depends on scaled parameters: the relative dielec- 
tric anisotropy a = Ae/e±, the Frank elastic constants 



Ki = K1/K3 and K2 = K2/K3 , 
and the magnetic field 

h = H/Hq with Ho = n2 — \ -r-^ . 

c y Ax 

If we introduce the magnetic coherence length 



6 



AxH^ 



(4.26) 



(4.27) 



(4.28) 



which gives the length scale over which director fluctu- 
ations are correlated, we obtain for the scaled magnetic 
field 



h = 



X 



27r6 



(4.29) 



where A = nii^jc. Thus, h ~ \ corresponds to a very 
short coherence length of X/2n. 

In the next subsection we will use the material pa- 
rameters of a typical nematic compound CBS to discuss 
the diffusion constants. 5CB is liquid crystalline at room 
temperature, and we use the parameters for 5 K below the 
nematic- isotropic transition [Q. The bend elastic con- 
stant is — 5.3 X 10^ dyne and is, as usual for conven- 
tional thermotropic nematic liquid crystals, larger then 
the splay and twist constants: Ki — 0.79 and if 2 = 0.43. 
For green light {lu/c = 1. 15x10'^ cm~^) the dielectric con- 
stant is e± = 2.381 and the anisotropy is a = 0.228 cor- 
responding to the refractive indices ni{d = 90°) = 1.710 
and — 1.543. Finally, the magnetic anisotropy is 
Ax = 0.95 X 10~^, from which we obtain a character- 
istic magnetic field of Hq — 3.6 x 10^ G. 



C. Light diffusion — results 

In this subsection we discuss the two essential compo- 
nents -Dm and D± of the diffusion tensor D{u}), which 



describe the diffusion of light, respectively, parallel and 
perpendicular to the director n p. B efo re w e can apply 
our general formulas from Eqs. (3.40), (3.41) and ( 3.49| ) 



we need to be more specific about our basis functions 
ipi(k), and we need the deriva tives of the components 



[Gq ]a which we gave in Eq. (2.2C). The calculations 



for the latter are straightforward and the results for the 
extraordinary and ordinary light mode are, respectively. 



ni(fc) 



and 



dk 



K 



k 

2 



C?cosi^fc 



CkK, 



(4.30) 



dk 



K 



2 



{\Jl- ClcoBLpkKi 



CkKu 



In the above, we chose 



(4.31) 



(4.32) 



wrote k in spherical coordinates as in Eq. (4.24) and then 
switched to the new C coordinate. The right-hand side 



of Eq. (4.31) for the ordinary mode is the isotropic result 
fc . K , which is modified in the extraordinary case. As 
our basis functions, we choose 



¥'7im(fe) = ipi^{Ck,ipk) 



(4.33) 



where ipimiC'k, ipk) stands for a real combination of the 
spherical harmonics Yim{Ck, ipk) and y/_m(Cfe, (^fc) in the 
new coordinate Ck instead of cosiJfc . It will soon become 
clear why the coordinate Ck is so helpful. The index 7 
stands for the basis of the tensor space, which we iden- 
tify here with our basis (g) 6^. Only functions with 
odd parity, i.e., with odd I, contribute to the diffusion. 
This is immediately obvious from the parity of the struc- 



ture factor and the derivative of \Gn 



Furthermore, 



we only need functions fimiCkj^k) containing cosmipk 
(to > 0); sinmipk is not necessary because of the sym- 
metry of [B^Q^^ (0)]q/3. With this choice of basis func- 
tions the matrix B is decomposed into submatrices for 



each TO, cx 6-^-^', because a term coarmpk cosmpk, 

(n ^ to), is not compatible with the rotational symmetry 
of the structure factor. Then, the diffusion constants _D|| 
and D± are, respectively, related to to = and to = 1: 



DnKi 



D±Ki 



with 



27l3 



2n^ 



J2[GiK)Uo[B-ryi'o[5iK)]yi,o (4.34) 



^ [g{K)],n [B'X-vi [^(^)]7'/'i (4.35) 



77'//' 
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n3 = nl{l + a/2) = (e±)^/'(l + a/2) . (4.36) 

1 are non-zero, 



Only elements of G{K) with I 
\Q{K )]yim 0, since the right-hand sides of Eqs. ( |4.30| ) 



and ( 4.31 ) depend only on I = 1 basis functions. This 
rises the question of how important higher I contributions 
are. For our calculations we, therefore, choose I — 1 and 



I = 3 functions 






fio{C, f) = 


VGttC 


(4.37a) 


ipiiiC, Lp) = 


VGtt \/i — coacp 


(4.37b) 


■^30 (C, (^) = 


y|7rC(5C2-3) 


(4.37c) 


(^31 (C, ip) = 


^^^Vl-C^l-5C^)cosip 


(4.37d) 



with the normalization 
dCdcp 



{2^f 



ipim{C, if) LpVm'{C,ip) =Su>S„ 



(4.38) 



The non-zero components of G{K) turn out to be 



[GiK)]no 
[GiK)]2io 
[G{K)]ui 



2e^{l + a)K\\/V6 
2ej_K\\/V6 

Kj^/Ve 



-2e_LvT 
[^(K)]2ii = -2e^KjV6 



(4.39a) 
(4.39b) 
(4.39c) 
(4.39d) 



The integrations to obtain the matrix elements of B re- 
quires more effort. We were able to calculate the integrals 
over if analytically with the help of 



dipkdfqfiifiq ~ tpk) = 2tt dipf{ip) 



dipkd(pqf{(p)cos^ipk = / dipf{ip) 



(4.40a) 



d(pkdipqf{ip)cosipk cosipq — dtf f [if)cofitf (4.40b) 



(4.40c) 



The results are listed in appendix C. The remaining in- 
tegrations were performed numerically. 

To discuss the diffusion constants we collect the pref- 
actors of the quantities involved to obtain an "averaged" 
transport mean-free-path 



c2 K3, 1 



(4.41) 



and write the diffusion constant Dn and ZJj^ in the form 



D\\=c^llD\\/'i , D^ = cxllD^/i 



(4.42) 



reminiscent of isotropic systems. The numerical factors 
Z?ll and D\_ only depend on the scaled parameters a, h, 
Ki, and K2, and c± is the speed of light of the ordi- 
nary light ray. We stress that /g is an averaged quantity 
and that our theory does not give a procedure to con- 
struct transport mean-free-paths for different light direc- 
tions (see however Ref. [^). The factor 97r in Iq is chosen 



such that Z?|| and D± are approximately 1 in the limit 
of an "isotropic" nematic with Ki — K2 = 1, a = 0, 
and /i = 0. We find b\\ = 1.053 and = 0.998 with 
a small anisotropy of D\\/D±_ = 1.06 because of the in- 
herent anisotropy in the nematic structure factor, which 
is re prese nted by the geometrical factor N{ea, e^, iis) of 
Eq. (4.21). We, at least qualitatively, understand why 
is larger than D± . The diffusion constants grow when the 
syste m's ab ility to sc atter light decreases. From equa- 



tions ( 4.34 ) and ( 4.35 ) for and D±, we recognize that 
the diffusion constants are, respectively, determined by 
scattering around the director (to = 0) or perpendicular 
to it (m = 1). However, forward and backward scatter- 
ing along the director is suppressed by the geometrical 
factor and we expect to be larger than D±. In a 
completely isotropic system, where we set N{ea,e0,us) 
equal to 1, the transport mean- free-path is easy to cal- 
culate. It is a factor of 4/9 smaller than Iq, which again 
demonstrates the effect of the geometrical factor. For 
temperature T = 300 K and the parameters of the com- 
pound 5CB, which we summarized in the last subsection, 
we obtain Iq — 2.3 mm in agreement with experiments 

In the following we will explore the dependence of the 
diffusion constants on the scaled parameters Ki, K2, a, 
and h. In Figs. || and ||we plot the relative changes of the 
diffusion constants when spherical harmonics of / = 3, in 
addition to ^ = 1, are included in the calculations. The 
field dependence in Fig. |5| shows that the changes are 
around 1 % or smaller and that D± is more strongly af- 
fected by higher spherical harmonics. The same is valid 
for the Ki and K2 dependence in Fig. |[ Only for ex- 
treme situations Ki < 0.1 or K2 < 0.1 do the changes 
grow to 3%. We conclude that the restriction to spheri- 
cal harmonics of Z = 1 gives a good approximation for the 
diffusion constants. The following graphs will, however, 
all be presented with the I — 3 contributions included. 

For the nematic compound 5CB we show in Fig. |^ 
how the diffusion constants Dy and D± and the rela- 
tive anisotropy (L)|| — D±)/D± behave in a magnetic 

field. and D± grow with H because the magnetic 
field suppresses director fluctuations. The field depen- 
dence of the relative anisotropy in the diffusion is weak. 
For ordinary magnetic fields up to 5 x 10^ G, which cor- 
responds to a magnetic coherence length ^3 of approxi- 
mately 1 fim, the changes in Z3j| and D± are small. The 
values ior H ~ read = 0.95 and D± = 0.65 with 
a ratio D\\/D±_ = 1.45 in good agreement with exper- 
iments The reason why we plot £)y, D± and 
(Z^ll — D±j/Dj_ on a large field range is to show that the 
quantities smoothly approach finite values at /i = 0. This 
is not obvious since the structure factor possesses a sin- 
gularity at /i = and for vanishing scattering vector q^. 
Strictly speaking, our weak-scattering approximation is 
not applicable here. However, in a completely isotropic 
model we understand why the quantities are finite. The 
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familiar factor 1 — cosz?^ 
from <7g cx 1 - 



in formula (3.53) cancels the 

COS!?,. 



singularit; 

In Fig. M we explore the anisotropy in the diffusion de- 
pending on the anisotropy a in the dielectric constants. 
As already discussed, when a = there is a small non- 
zero value of — D±)/D±. This grows with a because 
the speed of light of the extraordinary light ray is larger 
along the director than perpendicular to it. On the other 
hand, for a < —0.15, the anisotropy (-D|| — D±)/D± 
changes sign, and light diffuses faster perpendicular to 
the director. This effect and the inversion point -D|| = D± 
should be observable in discotic ncmatics where a is neg- 
ative. 




-0.012 



0.2 0.4 0.6 

i? /3.6 X lO^G 



0.8 



FIG. 5. Field dependence of the relative changes of diffu- 
sion constants after I = 3 spherical harmonics (superscript 
13) are included in addition to / = 1 spherical harmonics 
(superscript 1). Parameters are Ki — K2 = 1 and a = 0. 
The magnetic field is given relative to the characteristic field 
Ho = 3.6 X 10^ G of the nematic compound 5CB. 




FIG. 6. Relative difference as a function of K\ and K2 at 
a — 0, h = 0.1 between D*^' calculated with I = 1 spherical 
harmonics only and D^^' 
I = 3. 




0.4 0.6 0.: 

i? /3.6 X lO^G 



FIG. 7. Field dependence of the diffusion constants 73|| and 
D± and the relative anisotropy (Dy — D±)/D±_ for the ne- 
matic compound 5CB (Ki = 0.79, ¥2 = 0.43, a = 0.228 and 
Ho = 3.6 X 10^ G). 




calculated with both I 



1 and 



FIG. 8. Relative anisotropy (Dy — Dj_)/Dj_ depending on 
the dielectric anisotropy a for Ki — K2 ~ 1 and h — 0.01. 

Finally, we discuss the dependence of the diffusion on 
the elastic constants Ki and i^2- In Fig- H we show 
that D± decreases with the elastic constants since the 
light scattering from the director modes increases. At 
the extreme values Ki — K2 = 0.01, we find D± = 
0.07. The contour lines reveal an asymmetry between 
the splay (Ki) and the twist (^2) distortions. D± de- 
creases more strongly with K2. The diffusion constant 
_D|| shows a similar behavior. Fig. |l^ gives the anisotropy 
(£)|| —D±)/ D± for the same range. It grows with decreas- 
ing elastic constants showing that D± is more affected by 
splay and twist distortions than Dn . The asymmetry be- 
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tween splay and twist is clearly visible. The last two 
graphs cover the range of conventional thermotropic ne- 
matics where usually Ki < 1 and K2 < 1. In Fig. |ll| 
we extend this range to Ki = K2 = 10 and observe that 
the anisotropy {D\\ ~ D±)/D± changes sign. The contour 
line on the base of the coordinate system indicates where 
^11 = D±. Roughly speaking, < D± if Ki > 0.6 
and K2 > 1.4. In a smectic-A phase (SmA) twist and 
bend deformations are expelled by the layered structure 
pof , hence Ki <C 1. Unfortunately, this also means that 
certain scattering vectors show very weak scattering so 
that the diffusion approximation cannot be achieved for 
reasonably sized samples. However, in the vicinity of a 
SmA-nematic phase transition, where the layered struc- 
ture softens (SmA) or starts to form (nematic phase), the 
diffusion approximation of light could be used to study 
the behavior of the Frank elastic constants close to the 
transition. A second interesting system is a polymer ne- 
matic liquid crystal. For long rigid rods one expects a 
large splay constant ||6^,^. Taratuta et al. deter- 
mined the Frank elastic constants for a special system 
and found the ratios Ki = 0.85 and K2 — 0.07 with an 
absolute value of K3 — A.7 x 10"'' dynes, which is suit- 
able for the diffusion approximation of light. From these 
parameters we predict a "large" ratio of D^^/D± — 2.8. 
The reported system has a very low dielectric anisotropy 
a and scatters light only weakly. However, it should be 
possible to find systems which are more favorable regard- 
ing a. 

We also calculated a Taylor expansion for the relative 
anisotropy (_D|| — D±)/D± around Ki = K2 = 1 and 
a = and found 



^11 --Pi 



0.06 - 0.l{Ki - 1) - 0.3(^2 - 1) + 0.3a 

+ 0.1(7?i - 1)^ + 0.2(i?2 - 1)^ (4.43) 
-0.08(i?i - 1)(7?2 - 1) , 



where the coefficients are material independent. Second 
order terms in a and couplings to Ki and K2 are negli- 
gible. The expansion summarizes the whole discussion. 



D. Diffusing Wave Spectroscopy 

With Diffusing Wave Spectroscopy it is possible to 
measure the averaged dynamical properties of a system 
through the dynamic absorption coefficient /i(w, t) of Eq. 
(3.42). Direc tor modes are purely diffusive, as described 
by Eq. ( 4.20 ) for the structure factor, and possess a vis- 
cosity ris{qs) which we specify here |54|: 



{^i3ql - ^2gf )^ 



? 9 
^29^ 

V2{qs) = 7 2— 2 

Vaqi + Vbq\\ 



(4.44a) 
(4.44b) 




FIG. 9. Diffusion constant D± depending on Ki and K2- 
Parameters are a = and h = 0.1. 



Dii -D 



0.01 




0.01 



FIG. 10. Relative anisotropy (Dy — D±) / D± depending on 
Ki and K2, which range from 0.01 to 1. Parameters are a = 
and h = 0.1. 



Dii -D 




1 0.1 



FIG. 11. Relative anisotropy (Dy — D±) / D± depending on 
A'l and K2, which range from 0.1 to 10. Parameters are a = 
and h = 0.1. 
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with 



?yc = (a^s + M4 + M6)/2 

rim = fJ-l + Vb + Vc ■ 



(4.45a) 
(4.45b) 
(4.45c) 
(4.45d) 



The Leslie viscosities /Zi govern the viscous flow of the 
fluid and couple it to the director motion. The Miesow- 
icz viscosities ria, jyt,, and ri^ can be measured in pure 
flow experiments. The rotational viscosity 7 character- 
izes viscous forces due to rotations of the director. With 



the explicit formula (4.20) for the structure factor and 
small enough times to expand the exponential time fac- 
tor, ii{uj,t) becomes proportional to time 



/i(^,i)«(Ae)^fcBT^ ^ 



2"- Q,/3,5"''= 



miqs) 



t 



(4.46) 



If we collect all the prefactors, the dynamic absorption 
coefficient can be written as 



IJ,{uj,t) = Hot with Ho — a 



2 2k^T Lo'^ /i 



97r 



c\ 7 



(4.47) 



Here the numerical factor pi is a dimensionless angular 
average involving the geometrical factor and the viscosi- 
ties of the director modes. It depends on the Leslie 
viscosities relative to 7 and the dielectric anisotropy a, 
and it is equal 1 if ris{qs) = 7 and a = 0. When 
Vsils) — 1, H can be evaluated analytically even when 
a 7^ 0, and we find Ji = j^^yl- For the compound 5CB, 
7//I = 0.60 ± .20 poise was experimentally determined by 
Kao et al. using DWS and the last formula ||l|,||. This 
value is in good agreement with the rotational viscosity 
7 = 0.81 poise of 5CB and shows the validity of the 
theory. Of course 7 is larger than the measured 7//! since 
ris{qs) is smaller than 7 [see Eqs. (4.37)] so that Jl exceeds 
1. However, the values of the Leslie viscosities are such 
that Jl is of order 1 in usual thermotropic nematics. Fur- 
thermore, the Leslie-Erickson theory seems to describe 
the director modes properly down to 4 x 10~® s, the time 
resolution of the experiments, since there is no indica- 
tion for a deviation from the linear time depencence of 
/i(ci;,i) predicted by the diffusive director modes. Mate- 
rials where the viscosities ris{qs) for different values of 
Qs) differ by factors of 10^ — 10"^ are polymer nematic 
liquid crystals with some director modes relaxing on 
a much larger time scale than they do in ordinary ne- 
matics. It would be interesting to study such systems to 
see if they exhibit a deviation from the Leslie-Erickson 
theory for short times which, e.g., would show up in a 
different temporal power law for /x(w,t). 

Finally, we point out an important difference between 
nematic liquid crystals and coll oidal suspensions. The 
dynamic absorption coefficient (4.4f:) only contains the 



viscosities of the director modes. The Frank elastic con- 
stants cancel because they determine both the light scat- 
tering and the dynamics. On the other hand, in colloidal 
suspensions fiQ = 2Dbi^^ /{l*c), where the transport 
mean- free-path I* characterizes light propagation and the 
diffusion constant Db the Brownian motion of the col- 
loidal particles. 
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APPENDIX A: THE BETHE-SALPETER 
EQUATION 

The averaged two-particle Green function with all 
space-time coordinates looks like 

(G«(a;i,a=2;Ti,T2)®G-^(y2,yi;r2,ri))W . (Al) 

After an integration over electric field sources Wo(2) at 
points (a;2,T2) and (2/2,^^2) it provides us with the full 
autocorrelation function VF(1) of the electric light field 
at {xi,Ti) and (yi,Ti): 



W{1) 



rf2*(l,2)W^o(2) 



(A2) 



In the following we abbreviate a pair of points in space 
time by its index i. The two-particle Green function fol- 
lows from the Bethe-Salpeter equation 

*(1,4) = /(1,4) + I rf2d3/(l,2)t/(2,3)*(3,4) . (A3) 

The quantity /(1, 2) stands for the tensor product of the 
averaged one-particle Green functions: 

/(l,2) = [(G^)(a;i-a;2;ri-r2) 

®(G^)(y2-yi;T2-r)]W . (A4) 

It propagates two electric field modes from 2 to 1 between 
their scattering events. The irreducible vertex function U 
describes different characteristic sets of scattering events 
po| . We will soon specify it. Note, that there is no 
preferred point in time and that the scattering medium 
on average is homogeneous in space. Therefore, all our 
averaged quantities do not change under translations in 
space time, and hence they can only depend on differ- 
ences of the coordinates. We now introduce center-of- 
"mass" {Ri,Ti) and relative {ri,ti) coordinates: 
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Xi = + ri/2 

n = T, + u/2 



Vi = Ri 

f, = Ti 



-ri/2 
k/2 . 



(A5) 



It is straightforward to show that the Jacobian determi- 
nant for this coordinate transformation is 1, and we have 



di = d^R, d-^TidTdt 



(A6) 



All of our quantities only depend on differences in the 
center-of- "mass" coordinates Ri and Ti because of the 
homogeneity in space time. To discuss the Bethe- 
Salpeter equation we perform a Fourier transformation, 

/ d^{R^ - Rj) d^{T, - Tj) d^n d^rj d\ dHj 

. . . exp[-i{K ■ {R, - Rj) + k,-r,- kj ■ r/)] (A7) 
X exp[i(f2(Tj — Tj) + ujiti — iUjtj)] , 

which transforms our quantities as follows: 

*(1,4) ^ 



t/(2,3) 
/(I, 2) 

where = 



U' 



k2k: 



'{K,VL) 



(A8) 



= 27r5((jJi ~ ti'2), and Ife^^.^, defined in Eq. 
(3.29), contains 5{ki — ^2). The delta functions appear 
since /(I, 2) only depend s on the differences of the rela- 
tive coordinates. In Eq. ( |3.2^ ) we give the explicit form 
of f'^^{K,il). The Fourier transformed Bethe-Salpeter 
equation finally takes the form 



d^ki duji 
(27r)3 '2^ 



x*-,",(K,f]) = /-(K,^l)litU, 



(A9) 



and the autocorrelation function W{K, k, il, lo) for light 
with frequency oj and wave vector k follows from 



W{K,k,n,uj) 



d^k' dcj' 



(27r)3 2tt 



xWo{K,k',n,ij') . (AlO) 



So far, our manipulations are generally valid. Now 
we introduce some approximations. We use the weak- 
scattering approximation in which the irreducible vertex 
function is given by 

(7(1,2) «(<5£(a;i,Ti)®^£(yi,Ti))Wl|^^ 

xS{xi-X2)6{yi-y2)S{Ti-T2)S{Ti-T2) . (All) 

It only considers scattering events of the two electric field 
modes which are tied together through the structure fac- 
tor. Fig. n2 gives a graphic representation of U{1, 2). 



1 92 

Seixi, Ti) — ^;-^ S{ti - T2) S{xi - X2) 



U(1.2) = 



dr'i 



8{ti - T2)8{yi - 2/2) 



FIG. 12. The irreducible vertex function (7(1,2) in the 
weak-scattering approximation. The second pair of points in 
space time is tied to the structure factor via delta functions. 

In center-of- "mass" and relative coordinates we obtain 

J7(l, 2) = ( 5e(ri, ii) ® 5£(0, 0) ) W ^ ^ 

X 5{Ri - R2 + ^-^)6{Ri - R2 - '-^) (A12) 
X 5{Ti -T2 + i^)S{Ti - T2 - ii^) , 

where we neglect the partial derivative d/dT which 
probes time variations on much longer time scales than 
the time period of light. Then, in Fourier space, we ob- 
tain 



U^ltl^^K, r!) = -f {5e ® ^£)(^)(fei - fc2,c.i - C.2) , 

(A13) 

where the time derivative of a delta function was han- 
dled with its representation via plane waves. Since 
the temporal correlations in the dielectric tensor de- 
cay on a time scale much longer than the time pe- 
riod of light, U'^^'^^ {K , is strongly peaked around 
uji — L02. The Bethe-Salpeter equation (A9) then implies 
the same behavior for (/f , fi), i. e., #^^,'(K, fi) oc 
p{ijj — uj')g{ijj'), where p{lij — w') stands for the strongly 
peaked part around lj = to' and g{Lo') for the remaining 
smooth function in uu' . The Green function {K 
picks up a source term of frequency uj' and produces an 
autocorrelation function with frequencies lo centered nar- 
rowly around lo' . In the time domain this corresponds to 
exp(ici;'t) times a slowly varying factor in t. It is therefore 
appropriate to introduce the Green function 



d3(w-Lj')*fcfc''e-'("~"')*, (A14) 



which gives this factor for light sources of frequency uo' . 
If we multiply Eq. (A9) by exp[— i((jj — Lo')t\, rewrite the 
argument of {5e®5e)^'^ as uo—uoi = lo—oj'—{uji~oj'), and 
integrate over t we fi nally arrive at the Bethe-Salpeter 
equation of subsection IIIB for $^^,,(i<r, fi, i): 
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(2^) 



with 



ft'iK,il)lit (A15) 



^-^{Se®Se)^^\k^ki,t) . (A16) 



In deriving the last eq uatio n we replaced the argument 



Lo in f'^{K,Q,) in Eq. (|AS| ) by to' and the factor ujf in 

by (c^')^ 



APPENDIX B: TWO IDENTITIES 



1. Ward Identity 



than one. The last equation shows that the leading order 
in K and co mes so lely from Go{k,uj). 

In subsection III A we calculated the Green functions 
{G'^/^){k,uo) to first order in T,'^/^{k,Lo). They were 
diagonal and only involved the d iago nal elements of 
S^/'^(fc,a;). Let us look at Eq. ( |B5|) under this ap- 
proximation concentrating on the propagating part of our 
quantities only. As explained in the main text a Greek 
superscript or subscript corresponds, respectively, to the 
basis vector Baik) ® or d°'{k) (g) d". For /^(O, 0) we 
find 

ftiO, 0) = [(G^)(fe, <E> {G^){k, cj)] ^''^ 
«^[(G^)(fc,..)]"[(G^)(fc,c.)]^ 



a,f3 



ea{k)(g)e0{k)(g)ea{k)(g>ep{k) . (B6) 



The Ward identity establishes a linear relation between 
the irreducible vertex function U and the mass operator. 
In the weak-scattering approximation it says 

AS-(K,0) = J j^Bt,,{t = 0)AGt,{K,0) , (Bl) 

which can be proven by a variable transformation of the 
integrand. For the general case see VoUhardt and Wolfle 
ill. 



Terms like ei(X)e2(8iei(X)e2 appear because the suberscript 
(N) tells us to interchange the second and third basis 
vector in {G^) (g) (G^). They only couple non-diagonal 
components to each other. Since non-diagonal elements 
are beyond our approximation we don't have to consider 
them. Then we are able to write 



[fkioM 



a/3 



[(G^)(fc,c^r[(G^)(fc,c.r , « 

, a^P 



f3 



(B7) 



2. A useful identity 

In this subsection we derive a useful relation between 
AG%{K,n), f%{K,9), and AS^(K,rj). We start with 
the definition 

AGtiK,^) = {G"){k+,u;+) - (G^)(fe_,c._) (B2) 



The component refers to terms like ei (g) ei (g) 62 (S" 62. 
The derivation of Gq^ with respect to iv just gives 



(B8) 



The derivation with respect to k contains two contribu- 
tions: 



and insert the unit tensor 1 in the appropriate represen- 
tation: 

AGt{K,n) = (G^)(fc+,c.+)([(G^)(fc_,c._)]-i 

-[(G«)(fc+,c.+ )]-i)(G^)(fe_,c._) . (B3) 



Then we introduce (G^^^)(k+ uj+) from Eq. (|]TT]), 
use definitions ( |3.22D and ( p!2q ) for AS^(l^,f2) and 
f'^{K,Vt), and arrive at the final equation: 

AGt{K,n)=rk{K,n) [AT.t{K,n) 

-{Goi(fe+,c.+)-Goi(fc_,c._)}] . (B4) 

In the case ii', $7 0, we get to first order in K and Q.: 
AGt{K,n)^ftiO,0) [AS-(O,0) 

dGo\j^_dG^^^ ^ (B5) 



dk 



where we neglect first-order terms from and AS^ be- 
cause the components of AS^ are already much smaller 



dG^ 
dk 



J2[^^^^d^{k)^d^ik) 

a 

+[G„-i]„^K(fc)®d"(fc)] 



(B9) 



The second one only produces non-diagonal elements 
which we do not have to consider. This statement seems 
to be obvious because a small rotation of k rotates da(k). 
But our basis vectors are more general and we have 
to look at it more carefully. We have to show that 
ei{k) ■ ■£pd^{k) is zero. With = ep^i and Eq being 
symmetric and independent of k we can write 

e,(fc) • — d\k) = d\k) ■ — e,(fc) . (BIO) 
okj okj 

From the biorthogonality relation it is clear that 

d\k) ■ — e,(fc) = -e,{k) ■ — d\k) , (Bll) 
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which verifies the statement. We are now able to write 
down Eq. (|B5| ) within our approximation: 

[AG^(K,17r « [/r(O,0r" ([AS^(O,0)]. 



20 



APPENDIX C: STRUCTURE FACTOR 



We give the two important structure factors in scaled form and in the coordinates Cfe, Sk, Cq, Sq and (p. Prom the 
notation it is clear whether they belong to an extraordinary or ordinary light ray. The parameters are a, Ki, K2 and 
h. 



[-Bfei,2(0)]i2 



[-Bri,i(0)]ii 



1 



3. 

2 n2 



(l + a)2 Ql 



cos^t/? TVi + 2cos93 A^2 + 



K2QI + Ql+h^ 



with 



and 



<3n = {Ck - Cqf 



ATi = {SlCq - SlCk? 
N, = SlSl{Ck-Cqf 



The integration over ^p gives: 

j COSip [-6^1^2(0)] 12 dip ■■ 



2 SI 



1 1 



+ 



Q'i=Sl-2SkSqCosip + Sl 

N2 = SkSq{SlCq-SlCk){Ck-Cq) 
= {SlCq+SlCkf 

1 (I-^{Kr) I-\K2) 



K2 Ki C{h) \ Ki 



Kn 



+ 4SpiK2) 



77 5fe 

4 Si 



C{h) + 2K2{Sl - Si) C{h) + 2K^{Sl + Si) 



1 / C{h) + K,{Sl + Si) C{h) + K2{Sl + Si) 



C{h) 



K,I{K^) 



K^I{K2) 



_ _ ' 

+ ^[C{h) + K2{Sl + Sl)]I{K2) 



[Sfei,i(0)]iid^^ 



TT 

2 



C2 C2 

'-'k ,^2 , '-^f? ri2 
0„ + — Uj, 



c2 <i • c: 



1 1 

Wi~W~2 



2CkCq 



1 1 



+ 



C{h) 



if2 



+ ^ J(i^l) + SlCk) {Ck + Cq) + C{h)CkCq) 



J COSip [. 



Sfei,i(0)]nd^=- 



TT 1 



^ - SlCk)'C{h) + 2K,{StCl - S'qCl) {Si - SD) 



4 5353 

- i + ( C{h) + 2K2{Sl + S'q) 



+ 



1 ,.,2^ ^c:2^,2 /^£W+i^2(5|+5|) C{h)+K^{Sl + S', 



V ^2 ^(^2) 



K^I{K,) 



+ ASlSp{K,) 



C{h)+K,{Sl + Sl) 



(Cla) 
(Clb) 

(C2) 
(C3) 

(C4a) 



(C4b) 



(C4c) 



(C4d) 



(K^{SlCq + SlCk) {Ck + Cq) + C{h)CkCq) 



with 



C{h) = {Ck - Cqf + and I{Ki) = { ^,{Sk - Sqf + C{h)\ [^^(5^ + Sqf + C{h)\ }' 



-1/2 



(C5) 
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— UJ ^ ^ ~ UJ 

The structure factor [-B;.iqi (0)]ii diverges for /i = and q — > k . We give [-Bj.i i(0)]ii for small and 5C = 



(C6) 
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